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Abstract 

Criteria for the incipient motion of a rigid body initially resting on a rigid sur- 
face are formulated from first principles in this work. A modified Coulomb friction 
model and an associated distribution of reaction forces are proposed. There ex- 
ists a surprisingly large category of general motions, however, which subscribe to 
a more conventional analysis; an analysis made possible by identifying so-called 
"significant reaction surfaces". In this way a model which caters for the majority 
of combined translations and rotations is devised. Some introductry results demon- 
strate the accuracy with which fluids can be numerically modelled for the purposes 
of entrainment. This work is an extension of previous work by the same author. 
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There is nothing new about incipient motion analyses applied in a sedimentary context 
(see Buffington and Montgomery (1997)). Although the topic of sedimentation may not 
lend itself ideally to the application of an incipient motion analysis, it is an ideal context 
in which to formulate criteria for incipient motion, it being sufficiently general. 

A modified Coulomb friction model and an associated distribution of reaction forces are 
proposed in this work. There exists a surprisingly large category of general motions, 
however, which subscribe to a more conventional analysis; an analysis which depends on 
identifying so-called "significant reaction surfaces" . In this way a model which caters for 
the majority of combined translations and rotations is devised. 

The work presented in the subsequent sections is inextricably linked to the theory and 
model developed by the author in Childs (2000), Childs (1999) and Childs and Reddy 
(1999). This paper is a supplement and concludes part of that work. The numerical ap- 
proach overcomes at least two of the limitations faced by experimentalists. Localised 
intergranular flow is amenable to calculation (but not measurement). What was previ- 
ously an insurmountable problem in the laboratory, reduces here to an exercise in mesh 
generation. Shape can be more directly taken into account using numerical simulations 
and shape is important. 

Incipient motion assumes selective erosion is as, or more, important than selective depo- 
sition in an attempt to elucidate the problem of sedimentation. It denies the possibility 
that some deposits may exist in a state of dynamic equilibrium. 

2 Using Numerical Methods to Model Entrainment 
Forces 

The results for the "pebble in a pothole" (Fig. |I|) are intended to demonstrate something 
of the power with which state-of-the-art numerical methods are able to model fluid flow 
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and the forces exerted by it. The "pebble" (a die bead of neutral bouyancyQ) was released 
from rest at the centre of the standard driven cavity flow problem (see Fig. |2[) and its 
motion was accurately determined. 



The combined, dimensionless, free surface-fluid-rigid body problem can be stated as 
follows: Find x r b, 0, v r t>, uj, h and v (the dimensionless, respective position, orientation, 
velocity and angular velocity of the rigid body, the elevation of the free surface and the 
velocity field of the fluid), which satisfy 
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aF~ % J 



1 Although a neutral bouyancy is not immediately reminiscent of any real life sediment problem, it 

was used to give maximum meaning to the problem as a test. 

= I, J = 1 for correct implementations eg. backward difference for time integration - see Childs 

(2000) 
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Figure 1: The trajectories of various included rigid bodies released from rest at the centre 
of the driven cavity flow. Top Left: Re = 0.025, m = 251.3, J 33 = 314.2 and t = 3.6 
sees. Top Right: Re = 0.025, m = 251.3, J 33 = 1.0 and t = 4.0 sees. Bottom Left: 
Re = 0.025, m = 251.3, J 33 = 0.1 and t = 3.6 sees. Bottom Right: Re — 1, m — 1, 
moment of inertia (scaled) = 0.1 and t = 2.0 sees. 
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subject to the "no slip" requirements 

v \r rb = v rb + u A (x - c) 
at fluid-rigid body interfaces and 

v | r = 

at fixed, solid impermeable boundaries. Here Ju( no SU m) denotes the ith, dimensionless 
principal moment of inertia of the rigid body, H is the transition matrix for a transition 
to a reference whose axes coincide with these principal moments of inertia, c is the 
centre of mass of the rigid body, is the dimensionless surface of the rigid body, t is a 
dimensionless time, pf and p s are the fluid density and units of solid density respectively, 
m is the dimensionless mass of the rigid body, y%b is a dimensionless body force per unit 
mass, F is the deformation gradient, J its determinant, v mesh is the dimensionless velocity 
of a reference which is otherwise allowed to deform freely and Re is the Reynolds number. 
The Cauchy stress, which appears in the formula for P (the Piola-Kirchoff stress tensor 
of the first kind) is that given in terms of the constitutive relation for a Newtonian fluid, 
cr = —pi + 2pD. A fully viscous boundary layer is therefore an intrinsic property of this 
model. All quantities otherwise undefined have their usual meanings. 

The first four equations (Eqns. ([!]) - (|])) govern the motion of a rigid body in a fluid. 
Eqns. ([!]), (0) and (|3|) are a dimensionless form of Euler's equations written in terms 
of a flow fractional force. They embody conservation of angular momentum. Eqn. (^) 
embodies conservation of linear momentum. Eqns. ([5]) and (||) facilitate recovery of the 
solution ultimately sought, the orientation and position (as opposed to angular and linear 
velocity). Eqn. (^) is the free surface model, one in which surface tension is considered 
negligeable. Eqns. (2) and (9) are a completely general reference description of a viscous, 
incompressible fluid written in dimensionless form. It should, however, be noted that v is 
merely a function of a distorting reference and the description therefore remains inertial in 
the same way that Lagrangian descriptions do. Both fluid and free surface sub-problems 
are best solved using the finite element method. A fourth order Runge-Kutta-Fehlberg 
method was used to solve the entire, coupled rigid body sub-problem. Rigid body, free 
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Figure 2: The respective v± and i>2 profiles computed along two cuts through the centre 
of a driven cavity flow (see Childs (2000) and Childs and Reddy (1999)) for a range 
of Reynolds numbers up to 4000. 
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surface and fluid sub-problems were solved iteratively thereby enforcing the entire set of 
governing equations. 

A comprehensive derivation of all equations can be found in Childs (2000), Childs 
(1999) and Childs and Reddy (1999). Most of the numerical details, particularly those 
pertaining to the finite element method, can be found in the latter. Of course turbulence 
is not taken into account in this particular model. It is possible to incorporate turbulent 
features by way of modifying the viscosity according to a well established technique, the 
so-called k—e model. The model presented here may therefore be considered to have 
certain limitations. 

In the succesive Fig. [I] trajectories the mass was concentrated closer to the centre (a 
lower moment of inertia was used). The Reynolds number used was based on grain size 
in all but the bottom right set of results in Fig. [|. In this instance the so-called flow-field 
Reynolds number was used. 

3 Incipient Motion 

The possibility of mobilisation (or remobilisation) is, for sediments, perhaps more rel- 
evant than the transport of the rigid body by the fluid. Incipient motion can be used 
as a simplistic criterion on which to base deposition, consequently the hydrodynamic 
characterisation of sediments and their environments of deposition. 

3.1 Notation 

Suppose that the shape of the rigid body can be adequately described by some equation 
of the form 

f{x) = 



and that the rigid surface on which this rigid body rests can be adequately described 
in terms of some function b(xi,x 2 ) which specifies the vertical position of the surface in 
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Figure 3: T e and T c are the respective exposed and contact surfaces of a rigid body of 
shape f(x) = 0, resting on a rigid bottom of shape b(x\, x 2 ). 

terms of the horizontal coordinates X\ and x 2 . That is 

£ 3 = b(x 1 ,x 2 ) 

(see Fig. |3] for a schematic representation of the problem of interest). A contact surface 
and an exposed surface, T c and T e respectively, may be formally defined as 

r c = {x : f(x) = 0,x 3 = b(xi,x 2 ); x EH} 

and 

T e = {x : f(x) = 0, x 3 ^ b(x 1 ,x 2 ); x G TZ] 

in terms of this notation. The surface of the rigid body has outward and inward normals, 
n out and rij n respectively, 



I v/H 

One might expect f{x) and b(xi,x 2 ) to be complicated in real-life examples. 
3.2 Modelling Incipient Motion 

A more direct application of Newton's first law one possibly couldn't find than incipient 
motion. The force acting on a rigid sediment grain is best thought of as 
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where Ffl ow & g is the combined gravitational and flow-tractional force acting on the rigid 
body and F c is the sum of forces (reaction and frictional) acting at the bed-rigid body 
contact, T c . Any incipient translation will be initiated by a combination of gravitational 
and flow-tractional forces, 



modified by any induced forces of reaction and friction (cr is the Cauchy stress at r e , 
the rigid body-fluid interface, and mg is the weight of the particle). The combined 
gravitational and flow-tractional couple is 



The important thing to recognise here is that any incipient translation will nonetheless 
be a component of the force Ffl ow & g (other components contributing to couples or being 
cancelled outright) and any incipient rotation will likewise be a component of Tfl ow & g . 

A modified or extended Coulomb friction model can be arrived at by extending the 
model for rigid, planar contact surfaces to more general surfaces i.e. the distribution 
of the reaction forces along T c is obtained by similar use of Newton's third law. The 
reaction force for a plane contact surface is 




(10) 



' ' flow&g 





where c is some convenient pointQ and Q r b is the rigid body domain. 





and the associated frictional force is 




(13) 




and £ is the coefficient of static 



friction. Extending the above models to more general surfaces 
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2 The torque is independent of the point about which it was calculated 
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is obtained where H is the Heaviside step function and c\ is some scaling factor 

Cl Jr c H (n out ■ Ffl 0W & g ) dT c 

v ' 

C2 / the pertinent area 

in which c 2 is determined by requiring that a component of Ffl ow y g , normal to some part 
of the contact, be completely balanced by the same component of reaction force. The 
disturbing question of whether another such component will be likewise balanced, when 
using the same c 2 scale, then arises. This dilemma can be avoided for all but one category 
of incipient motion in the forthcoming analysis. 

For couples, 

T c = Ci [(x — c) A F 'reaction + {x — c) A Ff riction ] H (n out ■ Ffl ow y g ) dT c . 

There exists a surprisingly large category of general motions which subscribe to a more 
conventional analysis than the one just proposed. 

Possible modes of incipient motion are categorised as follows for the purposes of this 
work: 

1. Translation only (lifting or sliding). 

2. Sliding combined with an "away from the surface" rotation. 

3. Lifting combined with an "away from the surface" rotation. 

4. Pivotting combined with lifting. 

5. Rotation only (pivotting or an "away from the surface" rotation). 

6. Pivotting combined with sliding: 

(a) About the same contact. 

(b) About different contacts (less likely). 

By "away from the surface" rotations is meant that the sediment partcle is rotated off 
the contact surface in such a way that no frictional or reaction forces are incurred i.e. 
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Figure 4: The condition which amounts to immediately excluding axes of rotation around 
which normals to the contact surface can be drawn. 

there is no counter torque. This condition amounts (under all but the most exceptional 
circumstances) to immediately excluding axes of rotation surrounded and intersected by 
normals from the contact surface. In Fig. |] the axis about which an "away from the 
surface" rotation occurs lies in one of the unshaded regions - which one depends on the 
sign of the rotation. For flow-induced axes of rotation which lie in a region such as the 
shaded one in Fig. |], one might presume the point about which the rigid body will rotate, 
is the nearest point to the centre of mass which does not lie on a normal to the contact 
surface (based on a principle of least action). 



3.3 Incipient Translation Only {r flow&g = 0) 



steepest gradient 




Figure 5: The significant reaction surface to Ffl ow & g . 
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The significant reaction surface is that part of the contact surface which determines the 
direction of the incipient translation, here 

min n in ■ F flow&g . (14) 

To work out the nett or resultant force compute Ffl ow & g (Eqn. (fL0|)), then locate the 
significant reaction surfaceQin r c (Eqn. (H)). The resultant force which arises is 



F Fflow&g (F reaction F friction) H {ll>out ' F flow&g) (^^) 

where F 'reaction and F friction are given by equations ([12|) and flilf ) respectively. The re- 
sulting torque is 

T [(35 c) A F reaction {.X c) A -F/riction] H (^Tl ou f • F flow&g') 

where x* is the minimum prescribed by Eqn. fll4|). 



3.4 Incipient Rotation Only (F flow&g = 0) 



One can surmise that the torque Tfl ow & g is equivalent to a tangential force 

F T = T ^ gA(a? n7 C) . (16) 
1 1 a; — c 1 1 2 

acting at x, since Tfl ow & g = (x — c) AF T by definition. The counter-couple arising due to 
a reaction at the contact surface has a negative projection on the combined gravitational 
and flow-tractional torque at the desired location, furthermore, it is identified by the 
value of this projection being a minimum. That is, 



mm 

xer c 



f \ a ( T flow&g A (X - C) 

[x - c) A I - — K -r^ — - ■ rii 



'flow&g (17) 

A strategy for computing the torque exerted on the rigid body may therefore be sum- 
marised as follows. Compute Tfl ow & g (using Eqn. (|1~TD) then locate the position of the 
significant reaction surface in T c (using Eqn. (0))- The resultant torque which arises 
(assuming the minimum is unique for simplicity) can then be calculated using the formula 

t = T flow&g + {{n out ■ F T ) (x - c) A n in 

-C fn ction (x - c) A [F T - (n out ■ F T ) n 0Ut }}H (n out ■ F T ) (18) 



3 assumed unique for simplicity, otherwise the forces will be equally shared 
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steepest gradient 




Figure 6: The significant reaction surface to Tfl ow & g . 

where Cf r i Ct ion = rnin ( 1 , mt, — c -{ n out — \ Both frictional and reactional torques 

I \\* T-(n out -* -r)n out \\} 

are drawn up along similar lines to Ff riction and F reaction were. A non-zero value indicates 
the frictional force is insufficient to impede rotation. The translational force (a result of 
any uncoupling) is 

F = {{n out ■ F T ) n in — Cf riction [F T — (n out ■ F T ) n out }} H {n out ■ F T ) . 



3.5 Simultaneous Incipient Rotation and Translation (F flow&g ^ 

Assuming a common position for both the pivot and the initial sliding surface, it can be 
located using either the (|H) or flTT| ) minima. To determine the nett or resultant force 
and torque, compute Ffl ow & g and Tfl 0W & 9 (using equations fllO| ) and (|TlD ). Locate the 



position on T c at which either of the minima, (14) or (|T7|) , occur. At this point evaluate 
the quantities F T (using Eqn. QTB|) ) , 

Ffl ow g&T = Ffl ow &g + F T \X* , 



T^in ' F 'flowg&T \x* 



(19) 
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and 



[{x — c) A (Ffl owg& T ■ n in )] ■ Tfl ow y g \ x * 



(20) 



where x* is the prescribed minimum. The resultant force which arises in the event of 
Eqn. (|IgD being a negative quantity can be calculated according to 

F F 'flow&g {(flout ' F flowg&T) flin 

C 'friction [F flowg&T (flout ' F 'flowg &t) flout\} H (fl ou t ' F flowg &t) (21) 



where C 



friction 



min < 1 



C(n out -F 



out' 1 - flowg&T 



\F 'flowg&T- ( n ouf Ffl owg S-r) n out\ 



A zero value would corre- 



spond to no incipient motion. A non-zero value indicates the frictional force is insufficient 
to impede translation. The resulting torque is 

T — T flow&g + {{flout ■ Ffl owg&T ) (x* — c) A Tl in 

C friction [F 'flowg&T (.flout ' F flowg&T*) flout]} H (flout ' Ffl wg&-r) (^2) 

where x* is the prescribed minimum. The outcome could be summarised as follows: 



1. F j£ by Eqn. (0) and r = by Eqn. 

(a) minimum ( |T4"|) < => sliding 

(b) minimum ( |i~4"|) > =>- lifting. 



translation. 



2. Minimum (jl4|) < 0, minimum fll7|) > =>- sliding with an "away from the surface r 
rotation. 



3. Minimum (TA) > 0, minimum (|TT| ) > =>- lifting with an "away from the surface' : 
rotation. 



4. Minimum fllil) > 0, minimum ([[7]) < =r- pivotting. 



5. Minima (TTJ) and (JT7j) < 0, F = by Eqn. Q2Tj) and r ^ by Eqn. 
only. 



rotation 



(a) minimum (17) < =>- pivotting 

(b) minimum ( |T7j ) > =>- "away from the surface" rotation. 
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6. Minima (H) and (0)< 0, F ^ by Eqn. (0) and r ^ by Eqn. fl2|) slidin 
and pivotting. 

7. -F = by Eqn. (EHI) and r = by Eqn. (p2|) no incipient motion 



Note that "pivotting and sliding about different (less likely) contacts" is the only mode 
not comprehensively dealt with. 



4 Conclusions 

Classes of permissible incipient rotations and translations can be readily and systemati- 
cally formulated for rigid bodies placed on rigid surfaces, as can be a modified Coulomb 
friction model. These classes of permissible incipient motion together with the friction 
models, are proposed as criteria on which to base incipient translation and rotation. 

The possibility of mobilisation (or remobilisation) is, for sediments, perhaps more rel- 
evant than the transport of the rigid body by the fluid. Incipient motion can be used 
as a simplistic criterion on which to base deposition, consequently the hydrodynamic 
characterisation of sediments and their environments of deposition. The unknown final 
dynamic or static state of the depositional equilibrium could, in some instances, have 
serious ramifications for incipient motion based sedimentation models. 
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